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Abstract 

U 

We study the dynamics near finite-time singularities of flat isotropic universes filled with two 
interacting but otherwise arbitrary perfect fluids. The overall dynamical picture reveals a variety 
of asymptotic solutions valid locally around the spacetime singularity. We find the attractor of 
all solutions with standard decay, and for 'phantom' matter asymptotically at early times. We 
give a number of special asymptotic solutions describing universes collapsing to zero size and 
others ending at a big rip singularity. We also find a very complicated singularity corresponding 
to a logarithmic branch point that resembles a cyclic universe, and give an asymptotic local 
series representation of the general solution in the neighborhood of infinity. 



1 Introduction 

Fluid matter with all its ramifications has always played a key role in discussions of cosmological 
singularities. In studies of the genericity of quasi-isotropic solutions [1], in studies of the structure 
and nature of the singularity and energy conditions [2], in the construction of the general isotropic 
singularity |3] , in the singularity problem of inflationary cosmology , or in more recent attempts 
towards formulating the cosmological singularity in string and brane theory [5], one sees different 
manifestations of the 'nature abhors a vacuum' principle, i.e., using suitable 'fluids' to model the 
universe in its most extreme states. 
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In recent years there have been an increasing number of works devoted to analyzing diverse 
problems in situations involving more than one cosmological fluids that show a mutual interaction 
and the associated exchange of energy. Studies have been focused on a number of issues, for example 
a covariant description of the interaction [6], scaling solutions [7], perturbations [8], duality and 
symmetry transformations to obtain physically relevant solutions |9], detailed solutions with energy 
transfer |101 111), and of course on the important current issues of cosmic acceleration, dark matter 
and dark energy |12] . 

It is therefore important to understand the nature of finite-time singularities that may develop 
in cosmological models with interacting fluids. Such an understanding will complement current 
studies of such models which focus on other issues and may also provide a demarcation of the range 
of dynamical possibilities of these models. As this issue has not, to the best of our knowledge, been 
pursued in a systematic way so far, it is the purpose of this paper to carry out the first steps in 
providing the asymptotic properties of the solutions in the neighborhood of a finite-time singularity 
in cosmological models with two interacting fiuids. In particular, we shall focus exclusively on a flat 
FRW model containing two such fluids and construct asymptotic solutions which have the property 
to blow up at a finite-time singularity. 

The asymptotic analysis of the solutions of the dynamical system as the finite-time singularity 
is approached is carried out here using the method of asymptotic splittings, cf. [13^ 114] . In this 
method, the vector field that defines the system is asymptotically decomposed in such a way as to 
reveal its most important dominant features on approach to the singularity. This leads to a detailed 
construction of all possible local asymptotic solutions valid in the neighborhood of the finite-time 
singularity. These provide in turn a most accurate picture of all possible dominant features that the 
field possesses as it is driven to a blow up. For previous applications of this asymptotic technique 
to cosmological singularities, we refer to [T5]. 

The plan of this paper is as follows. In the next Section, we write the basic equations describing 
a flat FRW universe fllled with two interacting fiuids as a dynamical system and we are lead to the 
asymptotic fleld decompositions that will yield all possible dominant features as the singularity is 
approached. In the following Sections, we present an analysis of the asymptotic properties of the 
system generally divided into power-law, oscillatory and complete solutions. We conclude with a 
discussion in the last Section, pointing into more general aspects of this problem. 
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2 Asymptotic splittings 

We consider a flat FRW universe with scale factor a{t) containing two fluids with equations of state 

pi = iT-l)pi, p2 = i^-l)p2. (2.1) 

The unit and sign conventions we use are those of Weinberg [16], with SttG = 1. The total energy- 
momentum tensor is given by 

j.^(total)^^a)+T;2), (2.2) 

with Tfj:l) = (pi +Pi)u^Ui, +Pig^u-,i = 1,2, with = 5q, while we assume that the two fluids are 
not conserved separately, that is T^^^^^-u = ^ and T^'^^'^'^-y = —v^, so that j' (*°*^^) '^.'^ = g. 
The vector describes the energy transfer between the two fluids. The time component of the 
conservation equation, 

V^rW°'^=7;0, (2.3) 
(the spatial equations are trivial) becomes 

-a3pi + i-(a3(pi+pi)) = «V, (2.4) 
at 

or, 

^(pia3) + 3pia2 = aV, (2.5) 

or, finally, 

pi + ?,H{pi+pi) = Hav^, (2.6) 

with H = a/a for the Hubble expansion rate, and similarly for the second fluid. In general, we 
envisage an interaction of the form 

s = Hav^ = -I3H"'pI + aH^'p^, (2.7) 

where the exponents m, n. A, p are rational numbers indicating that between the two fluids there is 
an exchange of energy that depends nonlinearly on their densities and the Hubble rate. Depending 
on the signs of the constants a, /3, the fluids may 'decay' to each other transferring energy. Thus 
the evolution of this system is governed by the equations 

3H^ = PI + P2 
Pi + SHFpi = -/3H"'p^ + aH''p^ (2.8) 
P2 + 3HJP2 = pH'^p^-aH'^p^, 
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together of course with Eqs. (|2.1|) . (The first of these equations (Priedmann equation) is obtained 
as in the single fluid case, cf. [16j, p. 472, but with p = pi + P2 and p = pi+P2 in the expression 
for the total energy-momentum tensor.) 

Below we elaborate on the simplest case where the exponents m, n. A, p, are all set equal to onq^. 
This case corresponds to the problem studied in [10^ lllj . and it will be interesting to compare 
certain of our results with theirs. However, our approach is completely different, the focus here 
being exclusively on the asymptotic approach to the singularities of these models. Setting all 
the exponents m,n,X, p = 1 and renaming x = H, the system (|2.8p becomes equivalent to the 
dynamical system 

X = y (2.9) 

y = -Axy-Bx^, (2.10) 

where A = a + ^ + 3j + ST, B = 3{aT + + 3r-f)/2. This defines the vector field 

f{x,y) = iy,-Axy-Bx^). (2.11) 

The method of asymptotic splittings developed in [13^ [14] scrutinizes all possible modes that the 
vector field (j2.1ip attains on approach to the finite-time singularit}{^ located at t = 0. These modes 
correspond to the different ways that (j2.1ip splits as t — )■ 0. For the case we consider, the possible 
asymptotic modes are given by the following three distinct decompositions: 

f^^^ = {y, —Axy — Bx^), (all terms dominant case) (2-12) 
= {y,-Axy), (2.13) 
= iy,-Bx'). (2.14) 

Each one of the three decompositions ()2.12p - (j2.14p into which the vector field (j2.1ip splits, contains 
different dominant balances that describe the precise ways into which the dynamical system is 
driven asymptotically as we approach the singularity. These balances are in general non-unique. 
By further analyzing the balances of each particular decomposition, we are led to the construction 
of a number of possible asymptotic formal series valid locally in the neighborhood of the singularity, 
or in the neighborhood of infinity (the latter correspond to the behaviour of the system away from 



^The previously introduced more general couplings in Eq. (|2.7|) correspond to perturbations of this 'standard' 

case. It is not a priori obvious, however, that all such perturbations are physically relevant or distinct. 

^by a solution with a finite-time singularity we mean one where there is a time at which at least one of its 

components diverges. We note that the usual dynamical systems analysis through linearization etc is not relevant 

here, for in that one deals with equilibria, not singularities. 
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singularities, describing all possible complete solutions). From a close examination of the form of 
these asymptotic representations, we can obtain valuable information about the genericity of the 
asymptotic solutions, the stability/attractor properties of dominating solutions in the developments, 
and other precise information most valuable to create a detailed shape of the asymptotic evolution. 
We shall briefly comment on possible more general forms derived from the system (|2.8p in the last 
Section. 



3 Power-law solutions, the 6^0 attractor 

We start here our analysis of the possible asymptotic solutions towards the finite time singularity 
of the system (j2.8p by searching first for power-law type solutions. The first such solution we give 
in this Section is the simplest and perhaps the most important of them. Let us take the second 
decomposition 

/(2) = (y,-A:Ey), (3.1) 

and look for the possible dominant balances, by substituting in the system {x,y){t) = /(^^ the 
forms 

x{t) = etP, y{t)=Ct'', (3.2) 

where the coefficients 3 = {0, ^) G C, while the exponents p = (p, q) G Q. This leads to the unique 
balance 

fif) = [3,p] = [(2M,-2M),(-l,-2)], ^/O. (3.3) 
The candidate subdominant part /(^.^ub) _ —Bx^) of the vector field /^^^ satisfies 

(3.4) 



Here we have utilized the Barrow-Clifton parameter [10 



<5=^, (3.5) 

that will play an important role in the following. There is no way for the vector field /(^.^ub) j-^g 
subdominant asymptotically in the sense that 

^ ^p_i ' ^ 0, as t^O, (3.6) 

unless we set 

5 = 0. (3.7) 
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Otherwise the decomposition /^-^^ would not be acceptable asymptotically. This means that in 
order to satisfy this constraint, the subdominant part has to be vanishing. We take in this case the 
subdominant exponent q to be equal to one, cf. jl3] . 

Next we calculate the Kovalevskaya matrix (/C-matrix in short), given by 

/C = Z)/(2)(H)-diag(p), (3.8) 

where Df^'^\S) is the Jacobian matrix of /'■^•', at H, which in our case reads: 



= . (3.9) 




The next step is to calculate the /C-exponents for this balance. These exponents are the eigenvalues 
of the /C matrix and constitute its spectrum, spec(/C^^^). The arbitrary constants of any (particular 
or general) solution first appear in those terms in the asymptotic solution series whose coefficients 
Cfc have indices k = gs, where g is a non- negative /C-exponent. The number of non-negative /C- 
exponents equals therefore the number of arbitrary constants that appear in the series expansions. 
There is always the —1 exponent that corresponds to an arbitrary constant, the position of the 
singularity (here at t = for notational convenience). If the balance Bl is to correspond to a 
general solution, then it must possess a non-negative /C-exponent (the second arbitrary constant is 
the position of the singularity). Here we find 

spec(/cf)) = {-l,2}, (3.10) 

so that Bi indeed corresponds to a candidate general solution. Substituting the series expansions 

x = s°^oCiit^"S y = s~o92^^"^ (3.11) 

in the system ()2.8p and after some manipulations to determine the coefficients of the expansions 
recursively, we arrive at the following asymptotic solution around the singularity: 

X = jt-^ + C2it - ^cl^t"^ ■ ■ ■ , (3.12) 

while the y expansion is obtained from the above by differentiation. Note the arbitrary constant 
C21 appearing in this expansion signifying that this representation corresponds to a general solution 
(we need two for this, the second is the arbitrary position of the singularity). 

As a final test for admission of this solution, we use the Fredholm alternative to be satisfied 
by any admissible solution. This leads to the following compatibility condition for the positive 
eigenvalue 2 and an associated eigenvector, f2 = (1, 1): 

f J • f /C - ^/ j Cj = 0, (3.13) 



where I denotes the identity matrix, and we have to satisfy this at the j = 2 level. This gives 

C21 = C22, (3.14) 

and this is indeed true as found previously in the recursive calculation. It follows from Eq. (|3.12p 
that all solutions are dominated by the x = H ^ solution, that is the solution 

H^'^r^, or a(t)~t2M^ ^^3-^5^ 

is an attr actor of all smoothly evolving solutions at early times, assuming the weight-homogeneous 
Z*-^^ decomposition asymptotically. 

A comment about the results of [10] is in order. They find that at early times the attractor 
solution takes the form: 

aBc(t) ~ as t^O, (3.16) 

whereas we find the form ()3.15p . In terms of the parameter 5 defined in (13. 5p . their solution ()3.16p 
is given by 

aBc(t) ~ (i-'/l^l) , [0,1/8), (3.17) 

and we see that our solution (|3.15p includes the 5 = member of the one-parameter family of 5- 
solutions of the form (I3.17p . To enable the comparison, we note that when yl < in (j3.15p we have 
that our solution goes as t"^/'"^', and so the Barrow-Clifton family is asymptotic to our solution, 

aBc(i) ^ i"^^'^', as (5^0, (3.18) 

(in this case, the exponent of the a-Qc{t) solution in (I3.17P tends to 1). This result means that one 
branch of our solution (13.15P represents a limit function for the Barrow-Clifton family of (5-solutions 
(I3.17p . Since for the validity of the /^^^ decomposition asymptotically we were forced to take (5 = 0, 
we arrive at the interesting conclusion that the Barrow- Clifton solutions (j3.16p are all dominated 
by the solution (I3.15P in this casC. 

4 Phantom singularities 

Let us move on to the asymptotic analysis of the decomposition /^^^ = (y, —Bx^). There are two 

possible balances here but only one is of interest for the power-law solutions of this Section (we 

^We will comment later on the 5 = 1/8 limit of the 5-parametric family of solutions. For the moment we note that 
as it is expected from (13.171) . as 5 1/8 all these power-law solutions for small t will tend to zero (except of course 
for possible particular exact solutions, those with a smaller number of arbitrary constants than the general solution) 
and hence are expected to lose their significance asymptotically (this is like taking the limit limfc_^+oo ~ with 

ce (0,1)). 
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analyze the second balance together with other oscillatory solutions in Section 6). Substituting in 
the system (x,y)(t) = /^^^ the forms ()3.2p . we find that this balance is given by 

= [H,p] = [(±727^,^72/^), (-1,-2)] , B<0, (4.1) 

(the two branches give analogous results as we shall see). The candidate subdominant part 
y(3,sub) _ —Axy) of the vector field /^^^ satisfies 

y(3,sub)(-^^p^ 



(0,-^^0 = (0,0), (4.2) 



(3) 

i.e., it vanishes only when we set A = 0,9,(^ ^ 0. We note that the balance I3\ corresponds to the 
limit 

5 -00, (4.3) 

we expect that it refers to different parts of the (5-family of solutions than previousljQ. With the 
Kovalevskaya matrix of this balance having 

spec(/cf)) = (-l,4), (4.4) 

we find after some manipulation that the series expansion corresponding to this case is given by 
the form 

2; = ±y^i"' + C4it'TY^cliy— t^--- , (4.5) 

while the y expansion is obtained from the above by differentiation. 

We note here that although the dominant term in this expansion is the same as in (j3.12p . 
the whole formal expansion is a different one. The arbitrary constant C41 appearing in the series 
(j4.5p signifies that this representation corresponds to a general solution. Indeed, this becomes true 
since the compatibility condition for the positive eigenvalue 4 (with an associated eigenvector say, 
^2 = (1,3)), 

vj ■ (^IC - c, = 0, (4.6) 

at the J = 4 level gives 

3c4i = C42, (4.7) 

and this is true as it follows from the recursive calculation. 

It follows from Eq. ()4.5p that assuming the weight-homogeneous f^^^ decomposition asymptot- 
ically, all solutions dominated by the balance B\ (that is, those included in the family defined by 



^Indeed, this range of S means that the fluid parameters F, 7 cannot be positive simultaneously (this is shown in 
detail in Ref. [10]). hence the title of this Section. 
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(j4.5|) ) are attracted on approach to the singularity by the asymptotic solution x = H ^ -^t ^. 
That is, the dominating solution 

H^^r\ or a(t)~r2/iJ^ (^4 8) 

— i) 

is an attractor of all smoothly evolving 'phantom' solutions at early times. Other solutions, domi- 
nated by the second balance of this decomposition, are elucidated in the next Section. 



5 Decaying cosmologies and the borderline case 



We now focus on the asymptotic analysis of the all-terms-dominant case, that is the decomposition 
fW = —Axy — Bx^). The subdominant vector field is the zero field in this case, and there are 
two balances: 



B 



(1) 



B. 



(1) 



2B ' 2B 
' A-y/A^-^B -A + yjA^-^B^ 



2B 



2B 



,(-1,-2) 
,(-1,-2) 



(5.1) 
(5.2) 



Our analysis closely monitors the different values 5 may take and we focus in this Section exclusively 
on power law solutions, leaving the treatment of cyclic solutions for the next Section. Regarding 
the first balance of the /^^^ decomposition, the Kovalevskaya matrix is given by 



1 



1 



and we find 



-/i + 6 -^ + 2 



spec(/cS^^) 



where /i 



i + VT 



-1, 



(5.3) 



(5.4) 



As in this Section we restrict attention to power law asymptotic solutions, we examine the case 
6 = 1/8. Then we find 

spec(/cf)) = (-l,0), (5.5) 

with corresponding eigenvector 

vj = {l,-l). (5.6) 

The solution of the system is particular (only one arbitrary constant, cf. \13\ 114] for this terminol- 
ogy) and is given by 

x = H = ^t-\ or a~t^/^. (5.7) 
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We notice the two branches of this solution, one describing universes collapsing to zero size asymp- 
totically > 0), and the other ending at a big rip singularity (vl < 0). 

For the specific case 5 = 1/8, we note the general solution found in [17j is given by 



H"^ = a-^{c3 + C4lna). (5.8) 

If we set C4 = 0, then this is the same as the 1-parameter solution (|5.7|) found above. An exact 
solution identical to our solution (I5.7p was also found in |T0] . 

For the second balance of f^^^ decomposition, power-law solutions can be found for 6 = 1/8 as 
well as for the standard 'decaying fluid' range < < 1/8. The Kovalevskaya matrix is given by 



1 1 \ _________ , + 

+ 6 t + 2 



162^ = ( / ^ ^ I , where cj) = ^ — —, (5.9) 



with eigenvalues 



spec(4^))= (-1,^). (5.10) 



Further, we notice that the case (5 = 1/8 of this second balance has the same eigenvalues, eigenvec- 
tors and solution as the first balance of this decomposition. These are 



with corresponding eigenvector 



and solution 



spec(4'^) = (-l,0), (5.11) 



rij = (l,-l), (5.12) 



^t-\ (5.13) 



Let us now turn to the behaviour of the asymptotic solutions with standard decay, that is for 
< (5 < 1/8. For definiteness, we choose the value 5 = 1/9. Then 

spec(4'^) = (-l,l), (5.14) 

with corresponding eigenvector 

nj = (l,0). (5.15) 

After further manipulations, we find that in the asymptotic expansion arbitrary coefficients are 
expected to be in the places cn and C12, but from compatibility condition we have C12 = 0, giving 
therefore a solution with the correct number of arbitrary constants. The final solution is a general 
one and reads, 

3 A 

x = -r^ + cii + -cl^t + --- . (5.16) 
10 



The dominant part of the solution asymptotically is given by 

a ~ i^/^, (5.17) 
in accordance with the family found in |17j (for S = 1/9), that is 

= I + C^t + Cet^ (5.18) 
and it is also the member given in |10j for 5 = 1/9. 

6 Anti-decaying, cyclic and complete universes 

This section collects together all those cases where the asymptotic solution shows a qualitatively 
different character than that considered so far. The /^^^ decomposition gives imaginary solutions 
for B > 0. There are two balances: 

= [{iy^,-i^),{-l,-2)] (6.1) 
= [{-i^/2/B,i,/2/B),{-l,-2)]. (6.2) 

Upon considering the subdominant part, in the case i? > we find that this decomposition is 
asymptotically acceptable only if A = 0, therefore when 

S — )• oo. 

The eigenvalues of the Kovalevskaya matrix are for both balances given by 

spec(/cf^) = spec(/c5^^) = (-1,4), (6.3) 

with corresponding eigenvector 

uj = (l,3). (6.4) 

The coefficients 041,042 are expected to be arbitrary and the compatibility condition fixes one of 
them in terms of the other, 8041 = C42. The final solution is given by the expansion 

X = TiVVB t-^ + C4it^ T + ■■■ . (6.5) 

This solution has J — )• 00 and so it must belong to the family of antidecaying fluids considered in 
[10) . It is interesting that asymptotically the scale factor turns imaginary, perhaps an indication 
that the metric in this case becomes asymptotically Euclidean. In this case we may consider as 
'physical' that branch of the solution that S. W. Hawking would call 'compact', having zero size at 
the singularity, cf. [18]. 



11 



The /^^^ decomposition, on the other hand, leads to a very comphcated singularity for (5 > 1/8 
for both remaining balances. After setting 6 = 1/2, we get 

spec(/ci^^) =spec(/C^^^) = (-l,3^i\/3), (6.6) 

with corresponding eigenvectors 



{1,2^ iVs). (6.7) 



For both balances, the second eigenvalue of the /C-matrix has positive real part. The solution of 
the system then reads 

X = 1 , (6.8) 



or, in terms of the scale factor we find 



l±iV3 

ar^t^^. (6.9) 



This defines a multifunction on approach to the t = singularity which is obviously a logarithmic 
branch point admitting no Puiseux series representation. In this case, the scale factor never returns 
to its original value no matter how many times t loops around zero. 

Lastly, we give another sort of solution. When all the eigenvalues of the Kovalevskaya matrix are 
negative, the solution escapes away from the singularity towards infinity. For the /^^^ decomposition 
of the system, such a state appears when we examine the first balance with < 5 < 1/8. If we 
choose 6 = 1/9, then 

spec(4')) = (-l,-2), (6.10) 

with a corresponding eigenvector 

vJ = {l,-3). (6.11) 

To construct a suitable expansion for this case, we need to take the multiplicative inversqj to be 
equal to one, s = — 1, and the coefficients C21, C22 are arbitrary. Following the method of asymptotic 
splittings, we are led to the compatibility condition 3c2i = — C22, and finally the asymptotic solution 
valid in the neighborhood of infinity: 

X = jt-^ + ciit-'^ + C2it-^ + ■ ■ ■ . (6.12) 

This is a general solution valid away from any finite time singularity, showing a standard decay 
between the two fiuids. 



^that is the least common multiple of the set of subdominant exponents and the positive Kovalevskaya exponents, 

cf. [nnH]. 
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7 Discussion 



In this paper we provided a demarcation of the singular phenomena that emerge when we consider 
two interacting perfect fluids in a flat FRW universe. We have examined what happens when we 
take this system asymptotically to a finite-time singularity. We have found a number of regimes 
described by different asymptotic solutions - seven different behaviours in all. 

There is an asymptotic solution that acts as an attractor, a limit function to a wide family 
of solutions parametrized by the parameter 6. This solution is a member of a family of singular 
asymptotes that has the same number of arbitrary functions as the general solution, and attracts 
all these smoothly evolving solutions at early times in the 'direction' S ^ 0. There is an analogous 
behaviour for the so-called 'phantom' regime of asymptotic solutions. There are also decaying 
solutions collapsing to zero size, and decaying solutions to a big rip singularity, but these are of less 
generality than the afore-mentioned behaviour, valid for special values of S. The general solution 
towards the singularity with 'standard decay' (that is in the range 6 G (0, 1/8)) was also picked by 
our asymptotic method, and it was constructed for a concrete parameter value. We found solutions 
of the 'antidecaying' type that approach the finite-time singularity turning purely imaginary in the 
parameter limit S oo, these are perhaps more amenable to a quantum cosmological description. 
We also gave a very peculiar solution having a log-type branch point singularity describing a 'cyclic' 
universe at 'early' times. Lastly, we have given the behaviour of solutions away from singularities 
and towards infinity. 

The existence of the singular behaviours unraveled in this paper makes the dynamics of cos- 
mologies with two interacting fluids especially interesting on approach to their singularities, and 
the singularity in such models deserves to be further studied. One aspect of the problem that is 
currently under study is whether these forms of approach to the interacting fluid singularities are 
stable to perturbations of the m, n, A, ^ exponents away from the value one we considered in this 
work. This may demand a reformulation of the problem using more suitable variables. Another 
important issue that is also under examination is precisely how the inclusion of curvature alters 
the behaviours found in the flat case and whether new and distinct forms are possible. We plan to 
return to these more involved issues in the future. 
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